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ABSTRACT 

Here we discuss the properties of star cluster winds in the supercritical, catas- 
trophic cooling regime. We demonstrate that catastrophic cooling inhibits super- 
winds and after a rapid phase of accumulation of the ejected material within the 
star-forming volume a new stationary isothermal regime, supported by the ion- 
izing radiation from the central cluster, is established. The expected appearance 
of this core/halo supernebula in the visible line regime and possible late evolu- 
tionary tracks for super-star cluster winds, in the absence of ionizing radiation, 
are thoroughly discussed. 

Subject headings: clusters: winds - galaxies: starburst - methods: numerical 



1. Introduction 

Within the volume occupied by a star 
cluster (SC), the energy injected by stel- 
lar winds and supernova explosions is fully 
thermalized via random interactions. This 
generates the large central overpressure 
that continuously accelerates the ejected 
gas and eventually blows it out of the 
star cluster volume to compose a super- 



wind. In the adiabatic solution of Cheva- 
lier & Clegg (1985; hereafter referred to 
as CC85; see also Canto, et al. 2000 and 
Raga et al. 2001) temperature and den- 
sity present almost homogeneous values 
within the central volume, whereas the 
expansion velocity grows almost linearly 
from km s _1 at the center, to the sound 
speed (c) at the cluster radius r = Rsc- 
There is then a rapid evolution as mat- 
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ter streams away from the star cluster 
and the wind parameters (velocity, den- 
sity, temperature and pressure) soon ap- 
proach their asymptotic values: V w — >■ 
V Aoo = (2L 5C /M sc ) 1 /2 „ 2c(R sc ), Pw ~ 
r- 2 , T w ~ r- 4 / 3 and thus P w ~ r" 10 / 3 , 
where L^c* and M^c* are energy and mass 
deposition rates. 

Recent results on the outflows expected 
from SCs have led us to realize that the 
adiabatic steady wind solution proposed 
by CC85 and followers, becomes inappli- 
cable in the case of massive and concen- 
trated clusters. Radiative cooling strongly 
modifies first the temperature distribution 
predicted for adiabatic stationary winds 
(T w ~ r~ 4 / 3 ) bringing suddenly and within 
a small radius, the temperature down to 
10 4 K, restricting then the X-ray emissiv- 
ity of the winds to a volume much smaller 
than previously thought (see Silich et al. 
2003 and 2004; hereafter referred to as Pa- 
pers I and II). Also, as shown in Paper 
II, for more energetic clusters, strong ra- 
diative cooling promotes the sudden leak- 
age of thermal energy right within the star 
cluster volume itself, and for the cases 
in which the radiative losses there exceed 
30% of the stellar energy deposition rate, 
when cooling becomes catastrophic, then 
the stationary superwind solution is to- 
tally inhibited. In this latter case, the 
rapid drop in temperature within the SC 
volume, leads to a sudden drop in central 
pressure and particularly to a sudden drop 
in sound speed. This inhibits the fast ac- 
celeration predicted in the adiabatic solu- 
tion as required by the flow to reach its 
adiabatic terminal speed (vaoo ~ 10 3 km 
s~ 4 ). The sudden drop in sound speed up- 
sets also the balance, demanded by the sta- 



tionary solution, between the stellar mass 
input rate and the rate at which mat- 
ter can flow away from the cluster, ie: 

M sc = 4:irR 2 sc p sc (Rsc)c(Rsc), and this 
inevitably leads to mass accumulation. 

Here we show how the metallicity at- 
tained by the ejected matter from a coeval 
cluster, strongly affects the limits found in 
Paper II bringing the energy input rate fur- 
ther down to much smaller values and thus 
affecting even lower mass clusters. We also 
show how after the thermal stationary su- 
perwinds are inhibited, and after an in- 
evitably short phase of matter accumula- 
tion within the SC volume, a second sta- 
tionary solution can be found. This is only 
possible while the stellar UV flux ionizes 
the deposited matter to compose a dense 
supernebula with an isothermal (T ~ 10 4 
K), slow (t>oo ~ 50 km s" 1 ), stationary 
wind (see section 2). Finally, as the evo- 
lution continues and the stars producing 
the UV photon output evolve into super- 
novae, a third quasi-adiabatic stationary 
stage could become possible after a sec- 
ond phase of matter accumulation within 
the star cluster radius. This could in prin- 
ciple evolve into a massive, cold, neutral 
outflow, however, its maximum velocity is 
well below the escape speed and thus the 
ejected matter, unable to escape the grav- 
itational potential of the star cluster, will 
sooner or latter recolapse to cause a sec- 
ond major burst of star formation. Our 
conclusions and the observational proper- 
ties of the various evolutionary phases are 
given in section 3. 
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2. The properties of massive star 
clusters 

Superstar clusters recently found by the 
Hubble Space Telescope in a large variety 
of starburst galaxies (for a review see Ho 
1997 and the recent proceedings from "The 
formation and evolution of massive young 
star clusters", Eds. Lamers et al, 2004) 
present a typical half-light radius ~ 3 - 
10 pc, and masses that range from several 
times 10 4 M to several times 10 7 M (see 
Walcher et al. 2004). These are now be- 
lieved to be the unit of violent star forma- 
tion in starburst galaxies. 

The theoretical properties of massive 
bursts of star formation strongly depend 
on the assumed stellar evolutionary tracks. 
Synthesis models (e.g. Leitherer & Heck- 
man 1995, Leitherer et al. 1999) assume 
further a stellar IMF and an upper and 
lower mass limit for the star formation 
event. In this way one knows that a coeval 
burst of 10 6 M with a Salpeter IMF, and 
stars between 100 M and 1 M , would 
produce, through winds and supernovae 
(SNe), an almost constant mechanical en- 
ergy input rate of ~ 3 x 10 40 erg s" 1 for al- 
most 50 Myr, until the last star of 8 M ex- 
plodes as supernova. At the same time, the 
flux of ionizing radiation is to remain con- 
stant at ~ 10 53 photons s _1 for the first 3 
Myr, to then steadily fall off approximately 
as t -5 . The above implies that the HII re- 
gion phase will last for about 10 Myr, time 
during which the cluster photon flux will 
have decrease by more than two and a half 
orders of magnitude with respect to its ini- 
tial value. The above values scale linearly, 
throughout the evolution, with the mass of 
the assumed cluster and thus the HII re- 
gion phase is in all coeval cases four to five 



times shorter than the supernova phase. 

When dealing with the outflows gener- 
ated by star clusters, another important in- 
trinsic property is the metallicity of their 
ejected matter. This is a strongly varying 
function of time, bound by the yields from 
massive stars and their evolution time. 
Thus, once the cluster IMF and the stel- 
lar mass limits are defined, the resultant 
metallicity is an invariant curve, indepen- 
dent of the cluster mass. Here we con- 
sider coeval clusters with a Salpeter IMF, 
and stars between 100 M and 1 M , as 
well as the evolutionary tracks with rota- 
tion of Meynet & Maeder (2002) and an 
instantaneous mixing of the recently pro- 
cessed metals with the stellar envelopes 
of the progenitors (see Silich et al. 2001 
and Tenorio-Tagle et al. 2003, for an ex- 
plicit description of the calculations). This 
leads to metallicity values (using oxygen 
as tracer) that rapidly reaches 14 Z (see 
Figure 1), and although steadily decaying 
afterwards, the metallicity remains above 
solar values for a good deal of the evolu- 
tion (for more than 20 Myr), to then fall 
to the original metallicity of the parental 
cloud. One of the main effects of an en- 
hanced metallicity of the ejecta is to boost 
its radiative cooling and in such a case, 
massive clusters may inevitably enter into 
the catastrophic cooling regime, envisaged 
in Paper II, to then find their stationary 
superwinds totally inhibited. Relevant also 
is to note that at the end of the SN phase, 
the mass reinserted into the interstellar 
medium through winds and SNe amounts 
to about 40% of the original mass in stars 
(see Leitherer & Heckman 1995, their fig- 
ure 53). Here we address the fate of such 
a vast amount of high metallicity matter. 
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Fig. 1. — The metallicity of the matter 
ejected by coeval bursts of star formation. 
The metallicity (in solar units) of the mat- 
ter reinserted into the ISM by SCs is plot- 
ted as a function of the evolution time. 
The curve is derived from the metal yields 
of Meynet & Maeder (2002) stellar evo- 
lution models with stellar rotation, using 
oxygen as a tracer. The estimate assumes 
a Salpeter IMF and 100 M and 1 M 
upper and lower mass limit. The model 
also assumes an instantaneous mixing be- 
tween the newly processed metals and the 
envelopes of the progenitors. Note that un- 
der these assumptions the curve is indepen- 
dent of the cluster mass. 



mass clusters, appear far from the thresh- 
old line that separates the catastrophic 
cooling regime from the stationary super- 
winds, and thus are to generate station- 
ary quasi-adiabatic winds similar to those 
proposed by CC85. More energetic (or 
more massive clusters), as their energy in- 
put rate approaches the threshold line will 
produce strongly radiative stationary su- 
perwinds with considerably smaller X-ray 
emitting volumes. Finally, cases that lie 
above the threshold line in the catastrophic 
cooling zone, will radiate a large fraction of 
the energy input rate within the star clus- 
ter volume and thus will be unable to gen- 
erate a stationary superwind. 
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2.1. The limits between superwinds 
and supernebulae 

Figure 2 shows results from our self- 
consistent radiative code (see Paper II) in- 
dicating the limiting energy, as a function 
of the size of the clusters, at which radia- 
tive cooling becomes catastrophic for mat- 
ter with metallicities similar to those ex- 
pected for the gas emanating from massive 
stellar clusters. As pointed out in Paper II 
there are three different kinds of solutions: 
low energy input rates, produced by low- 



Fig. 2. — The threshold energy input rate. 
The limiting energy input rate above which 
catastrophic cooling inhibits the thermal 
stationary superwind solution for SCs with 
Vaoo = 1000 km s -1 , as a function of the 
size of the clusters and for two different 
abundances (solar and ten times solar) of 
the matter injected by the SCs. 

Figure 3 incorporates the run of the 
metallicity of the matter injected by winds 
and SNe as a function of time (Figure 1) 
to the results from our self-consistent ra- 
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diative solution (Figure 2) and plots the 
resultant location of the threshold line for 
clusters with a Rsc = 3 pc (solid line) and 
10 pc (dashed line), the typical range of 
sizes of superstar clusters, as a function of 
time. As mentioned above, clusters of any 
given mass will evolve at an almost con- 
stant mechanical luminosity, crossing the 
diagram from left to right (as the 10 6 M Q 
cluster shown in the figure as a solid ar- 
row). The figure confirms that low mass 
clusters will produce quasi-adiabatic sta- 
tionary winds. However, clusters with a 
mass in the range 3 xlO 5 M Q to 3 xlO 6 
M Q would cross twice the threshold line 
during their evolution and thus will have 
their thermal stationary winds fully inhib- 
ited for a good fraction of their evolution. 
More massive and compact clusters, those 
injecting more than say, 10 42 erg s -1 (see 
Figure 3), will be unable to develop the 
superwind outflows during their evolution. 

All of the latter clusters, unable to rid 
themselves from the continuous input of 
matter from their stellar sources, are to 
cause a rapid accumulation of the cold 
(T ~ 10 4 K) recombined ejecta, now sud- 
denly exposed to the UV photon flux from 
the most massive stars. As shown below, 
this leads to a dense and compact, high 
metallicity supernebula, able to establish 
an isothermal stationary HII region wind. 

2.2. The evolution of supernebulae 

For clusters that enter the catastrophic 
cooling regime, matter accumulation within 
the star cluster volume, would rapidly lead 
to larger densities, while the stellar UV 
photon flux independently will sustain the 
temperature at T ~ 10 4 by photoioniza- 
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Fig. 3. — The evolution of the threshold 
energy input rate. The continuous changes 
in the metallicity of the ejecta (Figure 1) 
shift the location of the threshold line. The 
figure shows the run of the threshold en- 
ergy input rate for massive clusters with 
3 pc (solid line) and 10 pc (dashed line) 
radii, respectively. The horizontal solid 
line across the diagram indicates the en- 
ergy deposition rate of a 10 6 M cluster. 
The scale for different mass clusters is indi- 
cated on the right-hand axis. The areas oc- 
cupied by the various stationary solutions 
described in this study are also indicated 
in the figure. 

tion. 

The density psc will increase until 
the stationary condition M sc = 47r_R| c 
Psc{Rsc)c{Rsc) (where now c(R S c) = 
Chii ~ 10 km s _1 ) is once again fulfilled. 
Given the drastic drop in temperature of 
the ejecta that results from thermaliza- 
tion followed by catastrophic cooling (say 
from 10 7 K to 10 4 K), the density psc will 
have to increase by one and a half orders 
of magnitude to compensate the drop in 
sound speed (~ y - 5 ) an d then meet the 
isothermal stationary condition. Mass ac- 
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cumulation will last for a short time (r) 

_ 47T ( P2 - Pl )R 3 SC VL 



2L S c 



10 V- (1) 



Once this happens, an isothermal, sta- 
tionary, photoionized wind will begin to 
emanate from the star cluster surface. 
The isothermal steady state flow equations 
within the star cluster radius (R < Rsc) 
are 



du w q m 1 + 3u 2 w /c 2 HII 



dr 3p w 1 - u 2 w /c 2 



HII 



Pw = 



3m,,, 



(2) 
(3) 



and outside of the star cluster (R > Rsc), 



du w 
dr 



Pw = 



1u w cfjjj/ll 



r 1 c2 HIl/ U w 



Anxu w r 2 ' 



(4) 



(5) 



where q m is the mass deposition rate per 
unit volume inside a star cluster, chii = 
(kT w / /ip) 1 / 2 is the isothermal sound speed, 
fi p = 14/23m H is the mean mass per par- 
ticle, k is the Boltzmann constant, and 7 
is the ratio of the specific heats. The wind 
central temperature is supported by pho- 
toionization at a T c ~ 10 4 K level, and the 
wind central density can be found by iter- 
ations from the condition that the isother- 
mal sonic point (u w = chii) acquires its 
proper position at the star cluster surface, 
Rsonic = Rsc (for more details see paper II). 
The isothermal wind temperature, density 
and velocity distributions are presented in 
Figure 4a-c (dashed lines). 

The resultant core/halo supernebula 
density distribution (see Figure 4) can 
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Fig. 4. — Supernebulae and massive slow 
neutral outflows. The stationary inner 
structure (temperature, density and veloc- 
ity, (a-c, respectively) as a function of R (in 
pc), adopted by the outflow from a 10 6 M 
compact (3 pc radius) cluster is shown for 
the supernebula (dashed lines) and massive 
slow neutral outflow phase (solid lines). 
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be approximated by a constant value n c 
within R,sc and by a power law (r~ 2 ) for 
r > Rsc- 111 this case, the number of re- 
combinations per unit time throughout the 
supernebula volume is: 

"Rsc 



dN r 



dt 



i-Rst 
\ 



n 



(r)r dr+ 



n 2 (r)r 2 dr 



16vr 



(3n 2 c R 



>3 

l sc 



1 - 



3R 



sc 



(6) 



where Rst is the Stromgren radius and (3 
is the recombination coefficient to all but 
the ground level. The Stromgren radius is 
then: 

3 Rsc 



R 



st 



4l 



167r/3n2_R| c 



(7) 



Equation (7) implies that the supernebula 
should be completely ionized by the star 
cluster UV radiation if the number of pho- 
tons per unit time exceeds the critical value 

N cn t = ^/3n 2 c R 3 sc . (8) 



The bright, centrally concentrated high 
metallicity nebula and its moderate wind, 
causing a stationary core/halo structure, 
are to be maintained for as long as they 
remain fully ionized. Note however that as 
the flux of UV photons from an instanta- 
neous, or coeval, burst of star formation 
declines drastically after ~ 3 Myr (see Lei- 
therer & Heckman 1995) it would eventu- 
ally fall below the critical value. The con- 
tinuously reduced number of ionizing pho- 
tons, unable to balance the large number 
of recombinations in the ionized volume, 
forces the ionization front to recede super- 
sonically towards the stars (see Beltram- 
etti et al. 1981). From equation (7) one 



can show that it will reach Rsc when the 
number of ionizing photons drops by a fac- 
tor of four below the critical value. At this 
time the central HII region will begin to 
loose its identity and uniform structure, to 
evolve finally into a collection of individ- 
ual ultracompact HII regions around the 
most massive members left within the co- 
eval star cluster. 

Taking into account that the superneb- 
ula central density p c = M sc / ' AhR 2 sc chii 
and that the star cluster mechanical lumi- 
nosity L sc = \M S cVl :00 , where V AiO0 is 
the adiabatic wind terminal speed in the 
hypothetical absence of radiative cooling, 
one can rewrite equation (8) in the form: 

4(3L 2 SC /r, \ 2 

^ ' crit 



ZnnlRscVX \chii 



pM 2 



sc 



(9) 



3nfi 2 a R sc c 2 HII 

where fi a is the mean mass per atom (fi a = 
14/1 Iitlh) and the supernebula isothermal 
sound speed chii ~ 11-6 km s -1 . In such 
a case, our standard 1O 6 M , R sc = 3 pc 
cluster and its supernebula will require of 
an N crit ~ 1.7 x 10 52 s -1 . This value is 
almost an order of magnitude below the 
maximum number of UV photons emitted 
initially by our example 10 6 M star clus- 
ter. The UV flux from the aging cluster 
would however reach the critical value af- 
ter 4 Myr and 5 Myr, depending if one 
assumes solar or 0.1 solar as the original 
cluster metallicity (see Leitherer & Heck- 
man, 1995; their Figure 37), and the UV 
radiation a couple of Myr after this will 
not suffice to sustain even the star cluster 
volume fully ionized. This brings the su- 
pernebula phase to an end, while restrict- 
ing the ionized volume to ultracompact HII 
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regions around the most massive sources 
left in the cluster. Without the sufficient 
UV photons to maintain the SC volume 
fully ionized and at a temperature ~ 10 4 
K, the stationary isothermal wind condi- 
tion is inhibited and a new phase of matter 
accumulation will start. 

Figure 5 displays the H Q luminosity of 
our 10 6 M cluster that injects 3 xlO 40 
erg s^ 1 within a Rsc = 3 pc and thus ini- 
tially while the metallicity remains below 
Z , it produces a strongly radiative sta- 
tionary superwind (follow the horizontal 
line in Figure 3). The stationary strongly 
radiative solution leads to a temperature 
~ 5 x 10 5 K at a distance of 17.5 pc 
and thus at this radius the streaming wind 
matter recombines and is exposed to the 
stellar UV photon output. Photoioniza- 
tion causes in such a case, a thin station- 
ary shell with the matter that continuously 
streams with large speeds (~ 10 3 km s _1 ) 
across the recombining radius and leads 
thus to a top-hat line profile (Rodriguez- 
Gonzalez et al. 2004; Owocki & Cohen, 
2001; Dessart & Owocki, 2002) shown by 
dotted line in Figure 5. This is very differ- 
ent to the calculated Gaussian line profile 
produced during the supernebular phase, 
which displays a linewidth of the order of 
100 km s _1 (FWZI) and, given the large 
densities, has an intensity several orders 
of magnitude larger than that of the top- 
hat line profile produced during the initial 
superwind stage. The total integrated in- 
tensity over the line profiles is: 5.3 xlO 40 
erg s _1 for the supernebula stage and 6.9 
xlO 34 erg s" 1 for the strongly radiative su- 
perwind. 



2.3. Gravitationally bound massive 
neutral outflows. 

The lack of sufficient UV photons to 
keep the supernebula phase at work, in- 
evitably leads to a further accumulation of 
the ejected matter within the star cluster 
volume, until the density there raises by 
another order of magnitude and compen- 
sates in this way the drop in temperature 
(say, from 10 4 K to 100 K) and a stationary 
flow can once again be established. This 
time however, the quasi-adiabatic flow is 
very massive and very slow (v^ ~ 2 km 
s _1 ~ 2cf/ 2 ). Figure 4 a-c (solid lines) dis- 
play the properties (the run of tempera- 
ture, density and velocity as a function of 
distance to the star cluster center) of such 
stationary flow. 

Figure 4 compares the isothermal (T ~ 
10 4 K) high metallicity supernebula out- 
flow with the late massive adiabatic neu- 
tral flow. In both cases the density drops 
as r~ 2 and thanks to the continuous en- 
ergy input rate through photoionization 
the supernebula is able to expand with in- 
creasingly larger velocities (up to 50 km 
s _1 ). The adiabatic dense outflow case 
reaches instead a maximum outflow veloc- 
ity of only v ~ 2c ~ 2 km s -1 . This maxi- 
mum speed is well below the escape speed 
(v esc = (2GM SC /Rsc) - 5 ) from the central 
cluster and thus is to evolve instead into 
another phase of matter accumulation, to 
eventually recolapse into a new generation 
of stars; subject of a forthcoming commu- 
nication. 
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3. Feedback and the observational 
properties of the superstar clus- 
ter outflows. 

Contrary to the predictions from the 
adiabatic solution of Chevalier & Clegg 
1985, we have shown here that massive 
and concentrated star clusters do not nec- 
essarily blow a strong stationary super- 
wind throughout their evolution. The high 
metallicity of the matter deposited through 
winds and supernovae, enhances radiative 
cooling within the star cluster volume, and 
leads instead to two other possible types 
of stationary flows emanating from mas- 
sive SC. These have been termed here su- 
pernebula and massive slow neutral out- 
flows. The high metallicity supernebulae 
with a distinct core/halo density distribu- 
tion, expanding with up to 50 km s _1 , can 
only be maintained in presence of an am- 
ple supply of stellar UV photons. These 
should be able to keep, at least, the central 
SC volume fully ionized. We have shown 
that this phase is thus bound to the first 
10 Myr of the cluster evolution. 

On the other hand, the slow massive 
neutral flows, that follow the supernebula 
phase, are to last in the case of very mas- 
sive clusters (M* > 3 x 10 6 M ), until 
the end of the supernova phase (~ 40 - 50 
Myr). The maximum speed of the neutral 
outflows is however much slower than the 
escape speed from the massive clusters and 
thus are to be inhibited by the potential of 
the clusters, leading instead to a further 
accumulation of matter which may even- 
tually collapse and lead to a new stellar 
generation. 

Very massive clusters, above the thresh- 
old line, will enter from the start of their 



evolution the supernebula phase, and will 
never cause a superwind. On the other 
hand, lower mass clusters will never expe- 
rience either the supernebula nor the mas- 
sive neutral stationary outflows and will 
establish stationary quasi-adiabatic super- 
winds throughout their evolution. Note 
however, that coeval clusters with a mass 
in the range 3 xlO 5 M to 3 xlO 6 M 
and a Rsc = 3 pc (see Figure 3) will start 
their evolution causing the development 
of a stationary strongly radiative super- 
wind, to then enter into the supernebulae 
phase and finally, after 8-20 Myr, respec- 
tively, re-enter once again into the station- 
ary strongly radiative superwind regime. 
Clearly, this will lead to a very strong in- 
teraction between the new superwind and 
the matter left near the cluster during the 
supernebula or the massive and slow neu- 
tral outflow phases. 

All superstar clusters will be, through- 
out their evolution, strong X-ray emitting 
sources. Clusters below the threshold line, 
that separates stationary superwinds and 
catastrophic cooling solutions (see Figure 
3), are to emit copiously in X-rays both 
within the SC central volume and in their 
winds. Clusters above the threshold line, 
whether experiencing the supernebula or 
the massive slow neutral outflow, will also 
emit in X-rays within their central volume, 
as catastrophic cooling radiates away the 
continuously replenished thermal energy. 
The X-ray properties of such clusters, com- 
pared to the X-ray emission from super- 
bubbles will be the subject of a forthcom- 
ing communication. 
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Fig. 5. — The signature of station- 
ary strongly radiative superwinds and su- 
pernebulae. The H a line profiles produced 
by the outflow from a compact (3 pc ra- 
dius), 10 6 M star cluster, at different 
stages of its evolution. The dotted line rep- 
resents the low intensity, flat-top H a pro- 
file calculated for a strongly radiative su- 
perwind stage, before the metallicity of the 
outflowing gas exceeds Z Q . In such a case 
strong radiative cooling brings the super- 
wind temperature down to T ~ 10 4 K at a 
distance of 17.5 pc allowing the stellar UV 
photon output to ionize the rapidly moving 
stream, while causing a broad flat-top line 
profile. This is to be compared with the 
line profile calculated for the supernebula 
phase (solid lines). The total integrated 
intensity over the line profiles is: 6.9 x 10 34 
erg s -1 and 5.3 x 10 40 erg s -1 , respectively. 
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